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ABSTRACT 



Aims. Our purpose is to place firm observational constraints on the three most widely used theoretical models for the spatial configu- 
ration of the large-scale interstellar magnetic field in the Galactic disk, namely, the ring, the axisymmetric and the bisymmetric field 
models. 

Methods. We use the rotation measures (RMs) of low-latitude Galactic pulsars and combine them with their dispersion measures and 
estimated distances to map out the line-of-sight component of the interstellar magnetic field in the near half of the Galactic disk. We 
then fit our map of the line-of-sight field to the three aforementioned theoretical field models and discuss the acceptability of each fit, 
in order to determine whether the considered field model is allowed by the pulsar data or not. 

Results. Strictly speaking, we find that all three field models are ruled out by the pulsar data. Furthermore, none of them appears to 
perform significantly better than the others. From this we conclude that the large-scale interstellar magnetic field in the Galactic disk 
has a more complex pattern than just circular, axisymmetric or bisymmetric. 

Key words. ISM: magnetic fields - Galaxy; disk - Galaxies; magnetic fields 



. 1. Introduction 

■ The interstellar magnetic field of our Galaxy has been the object 
of intense investigation since the early 1980s. Different obser- 
vational methods (e.g., based on synchrotron emission, Faraday 
rotation, Zeeman splitting, polarization of starlight, polarization 
of dust infrared emission) provide information on the magnetic 

■ field in different interstellar regions. Faraday rotation of Galactic 
' pulsars and extragalactic linearly polarized radio sources make 

it possible to directly trace the magnetic field in ionized regions. 
In practice, one measures the so-called rotation measure (RM), 

■ defined by 



directly yields the «e-weighted average value of B|j along its line 
of sight. 



^„ = 1-232 ^ .G 



(3) 



RM = 0.81 



B\\ ds rad m 



(1) 



We now know that the interstellar medium (ISM) is highly 
inhomogeneous and that the interstellar magnetic field has an 
important turbulent component. For this reason, neighboring 
pulsars may have significantly different values of RM and DM, 
and a plot RM versus DM will generally exhibit a large scatter 
However, if one considers a Galactic region larger than the scale 
of the turbulent field and containing enough pulsars for statisti- 
cal purposes, one can infer the large-scale (or regular) compo- 
nent of B\i in that region fro m the slope of the mean DM-RM 
relation dRand & Lvnelll994 : 



where is the free-electron density (in cm"^), By is the line- 
of-sight component of the magnetic field (in fiG) and d is the 
distance to the radio source (in pc). Pulsars present a number 
of advantages when used as probes of the interstellar magnetic 
field. In particular, they are highly linearly polarized, they have 
no intrinsic rotation measure and their distances can be estimated 
reasonably well. Moreover, the RM of a pulsar divided by its 
dispersion measure (DM), 



1-232 



dRM 
dDM 



(4) 
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Various theoretical models have been proposed to describe 
the spatial structure of the large-scale magnetic field in the 
Galaxy. First and foremost are the ring model, the axisym- 
metric or axisymmetric spiral (ASS) model, and the bisym- 
metric or bisymmetric spiral (BSS) model. According to the 
galactic dynamo theory, ASS fields would b e easiest to am- 
plify under typical galactic co nditions (e.g., iRuzmaikin et al.l 
ll985l:lFerriere & Schmitlll2000h . whereas BSS fields could pos- 
sibly be excited in the presence of an external disturbance, 
such as a companion galaxy (Moss 1995, 1996). On the other 
hand, the primordial field theory naturally leads to BSS fields 
(iHoward & Kulsrudlll997h . 
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In principle, RM studies are ideally suited to establish 
the overall structure of the Galactic magnetic field. However, 
the different RM studies performed so far yield contradic- 
tory results: some f avor a ring field (lRand& Kulkarni|[T98l 
iRand&Lv ni 11994 IValled l2005h . others an axisymmetric or 
ASS field (Valle Jll99ll Il996l). and others a bi s ymmetric or 
BSS field (Simard-Normandin & Kronberg '1980'; ' Han & Oiaol 
[1994; Indrani & Deshpande 1999; Han et al. 2006). Moreover, 
although all these studies conclude with a preferred field model, 
none of them has seriously considered the possibility that more 
than one model is allowed by the RM data or, alternatively, 
that none of the three basic models alone can account for the 
data. Hence the question we would like to address in this paper: 
which among the ring, axisymmetric and bisymmetric models 
can clearly be accepted on the grounds that it is consistent with 
the RM data, and which model should clearly be rejected on the 
grounds that it fails to provide a good fit to the data. 

In recent years, numerous pulsars were discovered in the 
near half of the Galactic disk and many of them had their RM 
measured. At the present time, among the ~ 1800 known pul- 
sars, 690 have measured RMs and, among the latter, 524 are lo- 
cated at low Galactic latitudes (\b\ < 10°). Pulsars with measured 
RMs now provide a reasonably good coverage of the near half 
of the Galactic disk. Furthermore, pulsar distances can now be 
estimated with fairly good accuracy thanks to the improved free- 
electron density model of Cordes & Lazio (2002) (known as the 
NE2001 model). The new measurements enable one to investi- 
gate the configuration of the Galactic magnetic field over a much 
larger region and with much more confidence than previously 
feasible. 

In Sect. 2, we present the three basic theoretical models for 
the interstellar magnetic field in the Galactic disk. In Sect. 3, we 
describe the procedure used to bin the pulsar data and to map out 
the distribution of By. In Sect. 4, we fit our map of By to each of 
the three field models, and we discuss how good the fits are at 
reproducing the pulsar data. In Sect. 5, we summarize our results 
and conclude our study. 

2. Description of the field models 

Throughout this paper, the Galactocentric cylindrical coordi- 
nates are denoted by (r, 0, z), and the distance from the Galactic 
center (GC) to the Sun is set to = 8.5 kpc. 

In general, the horizontal position of a given pulsar P can be 
defined either by its distance from the GC, r, and its Galactic 
azimuthal angle 6 (which increases clockwise from 6-0 along 
the line segment GC-Sun), or by its distance from the Sun, d, and 
its Galactic longitude, / (which increases counterclockwise from 
/ = along the line segment Sun-GC). Another useful angular 
coordinate is the angle a between the azimuthal direction at P 
and the vector P-Sun, such that a - 6 + I + ^ (see Fig.[T]l. 

Because Galactic differential rotation efficiently stretches 
magnetic field lines in the azimuthal direction, {Bg} dominates 
over both (B,) and (B,). Moreover, all the pulsars selected for 
the present work lie at \b\ < 10° and reside in the Galactic 
disk. There, the large- scale magnetic field is nearly horizontal 
(e.g.. lRuzmaikin et al.|[l985; Beck et al.,1996) , so that | (B^) | <k 
I (Br) I, I (Bg) I . In addition, projecting an already small (Bj) onto 
the line of sight to a pulsar further reduces its contribution by 
a factor \smb\ «; 1 . Under these conditions, the line-of-sight 
component of the large-scale magnetic field depends only on its 
radial and azimuthal components, and is related to them through 



We now present the three theoretical field models. 



2.1. Ring model 

In the ring model, the large-scale magnetic field points every- 
where in the azimuthal direction, so that its radial component 
vanishes: 



(Br) = . 



(6) 



Its azimuthal component is constant along circles, i.e., indepen- 
dent of 0, but it can vary with r : 



{Bg}^{Bg}(r), 



(7) 



and it can even change sign along a Galactic radius. As a matter 
of fact, all R M studies leading to a ring model have found rever- 
sals in (Bg) (iRand & Kulkarnilll989tl^nd & Lvndl 19941 1 Vailed 
2005). It should be noted that the ring model constitutes a par- 
ticular case of the axisymmetric model. 

2.2. Axisymmetric model 

In the axisymmetric model, (B,) and (Bg) are both independent 
of and vary only with r : 

{Br)^{Br)(r), (8) 



{Bg}^{Bg}(r) 



(9) 



Here, too, (Bg) can reverse sign with r. Such sign reversa ls were 
found in RM studies favoring an ASS magnetic field (IValled 
119911 119961) . Interestingly, reversals in (Bg) are also consistent 
with dynamo theory, which can produce them under certain 
conditions, e.g., when the magneto-ionic disk has a particular 
shape and thicknes s and the seed field itself has strong reversals 
(iPoezd et al.ll 19931) or when the Galactic rot ation rate decreases 
not o nly with radius but also with height jFerriere & Schmiti 
I2000h . 



2.3. Bisymmetric model 

In the bisymmetric model, (B^) and (Bg) have a simple si- 
nusoidal dependence on 6, which can be written in the form 
(iBerkhuijsen et al.lll997h : 



{Br)^b,.(r) sin(0-0(r)) 



(Bg) = bgir) sm(6 - cf>(r)) , 



(10) 



(11) 



where br(r) and bg{r) are the maximum amplitudes of (B, ) and 
(Bg), respectively, and <p{r) is the azimuthal phase. Both (B, ) and 
(Bg) can reverse sign with r. The magnetic pitch angle is defined 

as 



tan p(r) 



(Br) _ bAr) 
(Bg) bg{r) 



(12) 



(By) = (Br) sin a + (Bg) cos a , 



(5) 



it is positive (negative) if the magnetic field spirals out clockwise 
(counterclockwise) or spirals in counterclockwise (clockwise). 
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Fig. 1. Schematics showing the geometrical variables associated with a pulsar P (see main text for the exact definitions). 



3. Mapping of 

To date, there a re 69 
jHamilton & L vnd 11987 




pulsars wi th mea s ured RMs 



1997; Han etal. 



Weisberg et al 



Rand & Lvnd Il994t lOiaoetal 



..1999: iCrawford etal 
2004tlHan et al.ll2006t 



JMoutsos et al. 2008|). Among these pulsars, we selected those 
that lie at low Galactic latitudes {^\ < 10°) and have reliable 
RMs (error on RM < 25 rad m"^). This left us with 482 pulsars. 
For the distances and DMs of our selected pulsars, we used 
the va lues given in the ATNF Pulsar Catalog (Manchester et al.' 
|2005l see http://www.atnf.csiro.au/research/pulsar/psrcat). 
Pulsar distances in this catalog were estimated with the help of 
Cordes & Lazio's (2002) NE2001 model for the free-electron 
density; for pulsars located in the inner Galaxy, individual 
distances are typically uncertain by ~ 20%, but the relative 
distances of neighboring pulsars have a much lower uncertainty. 
Pulsar DMs, for their part, are known with good accuracy (error 
on DM generally < 1 cm"^ pc). 

In order to map out the large-scale component of By, one 
needs to divide the Galactic disk into regions (boxes) having 
sizes intermediate between the large scales of the regular field 
and the small scales of the turbulent fiel d and containing at 
least a few pulsars each. In previ ous studies jRand & Lvnell994t 
IWeisberg et all 120041: iHan et alJ[2006) . this division was based 
on a heliocentric grid defined by circles of constant d and ra- 
dial lines of constant I. Such a heliocentric division was jus- 
tified by the spatial distribution of the available pulsars, but 
it is ill-suited to the present work, whose purpose is to test 
field models expressed in terms of Galactic radius, r. A much 
more appropriate division here is one based on a hybrid grid 
defined by circles of constant r and lines of constant I (see 
Fig. |2]i. To make full use of the pulsar data, we consider two 
different grids. In the first grid, the circles are located at r = 
4 kpc, 5 kpc, 6 kpc, 7 kpc and 8 kpc [i.e., r - r,, with 
r, = / kpc, / = 4... 8], and the lines of constant I are the 
lines emanating from the Sun and tangent to one of the cir- 
cles r - ri, i - 2...1, plus the line Sun-GC [i.e., I - k, l±i, 
with /o = and l-tt = +asin^, / - 2. ..7] (see Fig. |2^). The 
second grid is defined in an analogous manner with the circles 
shifted by 0.5 kpc [i.e., r - r,-, / - 3.5, 4.5 ...8.5, and accord- 
ingly, / = Iq, l±i, i = 1.5, 2.5 ...7.5] (see Fig.|2]3). To ensure a 
sufficient number of pulsars per box, some of the boxes defined 
by these grids are paired together More specifically, the non- 
outermost boxes along each ring are paired either with their left 
or right neighbor along the same ring (thereby leading to a sin- 
gle double-size box) or with both neighbors separately (thereby 
leading to two overlapping boxes). 



Altogether, we have 9 rings, which are centered on r, = i kpc, 
with i = 4.5, 5.5 ... 7.5 in the first grid and / = 4, 5 ... 8 in the sec- 
ond grid. In the following, the ring centered on r, (and extending 
between r,_o.5 and r,+o.5) is referred to as ring /. 

We retain only the boxes containing at least 5 RMs. This 
minimum number of RMs, which is smaller than generally rec- 
ommended for statistical testing, results from the limited num- 
ber of pulsars with know RMs; it was chosen as a trade-off be- 
tween the need to have enough data points per box to make use 
of Eq. (|4]i and the need to have enough boxes on the Galactic 
plane to capture the spatial variations of the large-scale mag- 
netic field. Even with such a small number of RMs in a given 
box, it is possible to trace the dependences of R M on distance 
and on DM for the pulsars lying in this bo x (see iRand & Lvnd 
n^lWeisberg et al.ll2004tlHan et al.ll2006h . 

We estimate the average value of By, denoted by (fiy), in 
each of the retained boxes in the following way: we plot the 
points (DM,RM) of all the pulsars in the box, we fit a straight 
line through the resulting set of data points, and we take (^B\\j to 
be 1.232 times the slope of this line (see Eq. (|4|i and preceding 
comment). To perform the straight-line fit, we resort to a slightly 
modified version of the ordinary le ast-squares linear regression 
of F on X described bv llsobe etall (|l990), which is well suited 
when the dispersion of the data points about the linear relation 
cannot be calculated beforehand. The modifications brought to 
the original scheme are designed to exclude t he occasional out - 
liers - such as those arising from Hii regions (Mitra et alj|2003h . 
In practice, we discard all the data points whose absolute devi- 
ation from the mean RM or DM exceeds three times the mean 
absolute deviation. To illustrate the procedure, we show two ex- 
amples in Fig. [3] 

The derived values of (^Biij in all the boxes of our two grids 
are mapped in Fig. H] For convenience, these values are con- 
verted into vectors oriented along the local line of sight. 

As explained above Eq. (HI, the turbulent component of the 
magnetic field, SB, causes the RMs to scatter about the mean 
DM-RM line. This physical scatter due to turbulence is typically 
one order of magnitude larger than the observational scatter due 
to measurement errors. The amplitude of the RM scatter is given 
by the r.m.s. deviation of the measured RMs from the mean DM- 
RM line. Since the RM scatter is of predominantly turbulent ori- 
gin, its amplitude divided by the mean DM in the considered 
box directly yields (to a factor 1.232) an estimate for the r.m.s. 
value of the turbulent component of the line-of-sight field, 68^. 
Finally, the rm.s. value of 6B\\ divided by the square root of the 
number of pulsars in the box provides an estimate for the sta- 
tistical uncertainty in cry, which, again, is predominantly 
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X (kpc) 

Fig. 2. Grids used in our study overlaid on the face-on distribution of our 482 low-latitude (\b\ < 10°) pulsars. Pulsars with a 
positive (negative) RM are denoted with crosses (circles), (ji:,^) are Galactocentric cartesian coordinates. The Sun is located at 
(x = 8.5 kpc,y - 0). In the upper panel, the circles underlying the grids are at r = 4 kpc, 5 kpc, 6 kpc, 7 kpc and 8 kpc, while in 
the lower panel, they are at r = 3.5 kpc, 4.5 kpc, 5.5 kpc, 6.5 kpc, 7.5 kpc and 8.5 kpc. 



due to turbulence . The exact expression of cr|| can be found in 
llsobe et all (Il990l) . With our data, the typical values of cr|| lie 
between ~ 0.2 fiG and 1 .3 /iG. 



4. Data fitting to the field models 

Once we have obtained a set of observational values of ^By^ to- 
gether with their statistical uncertainties (or error bars), we can 
put the three theoretical field models presented in Sect. 2 to the 
test. As in all other studies based on RMs, we proceed on the no- 
tion that the large-scale interstellar magnetic field may be iden- 
tified with its «£ -weighted average value (denoted with an over- 
bar). Implicit here is the assumption that fluctuations in mag- 
netic field strength and in free-electron density are statistically 
uncorrected. In reality, t his assumption is certainly not strictly 
satisfied in the ISM (e.g.. lBeck et al.ll2003h . and this will cause 
our results to be somewhat biased. 

With this caveat in mind, we now describe the overall pro- 
cedure. For each model, we use all our observational values of 
to derive the best-fit parameters of the model. We then ex- 



amine whether the best fit is consistent with the pulsar data, i.e., 
whether the theoretical line-of-sight fields predicted by it fall 
within the error bars of the observational ^By^ (in a statistical 
sense). If we find that the best fit is not consistent with the data, 
we may conclude that the considered model must be rejected. If, 
on the other hand, the best fit is found consistent with the data, 
we may conclude that the model is acceptable; we then deter- 
mine the extent of the so-called "consistency domain", i.e., the 
parameter domain around the best fit within which solutions are 
consistent with the data. 

It is important to realize that the concept of acceptability dif- 
fers from the concept of detectability. A given field model, say, 
the ring model, is acceptable only if it is not ruled out by the 
available pulsar data. This does not necessarily imply that the 
Galactic magnetic field is really of the ring type, nor that a ring 
field has truly been detected. Detection of a ring field requires 
not only that the ring model be acceptable, but also that the zero- 
field solution do not belong to the consistency domain. 

Let us now discuss more specifically what exact criterion 
should be used to test consistency with the pulsar data for a given 
field model. Each of the three models is characterized by a num- 
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X (kpc) 

Fig. 4. Face-on map showing the average Hne-of-sight component of the magnetic field, obtained in the different boxes of 

our two grids. Each ^By^ is plotted in the form of a vector centered on the box midpoint (point at middle radius, r - ri, and middle 
longitude, I - (/min + /max)/2), and oriented along the local line of sight. 
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Fig. 3. Plots of RM versus distance (left panels) and versus DM 
(right panels) for the pulsars lying in two different boxes. The 
first box (top row) is delimited by the circles r - rg 5 and r = re.s 
in the radial direction and by the tangential lines I - /g 5 and 
I - 4.5 in the longitudinal direction. The second box (bottom 
row) is delimited by the circles r = rg 5 and r = r7.5 and by the 
tangential lines I - U,^ and / = /e.g. For each box, the best-fit 
straight line through the points (DM,RM) is drawn in the right 
panel, and the corresponding value of (^B\\^ with its statistical 
uncertainty are written in the upper right corner. The red points 
represent outliers. 



ber of independent free functions of Galactic radius [{Bg} (r) in 
the ring model; {8^} (r) and (Bg) (r) in the axisymmetric model; 
br{r), bg(r) and (p(r) in the bisymmetric model], corresponding 
to the same number of independent free parameters in every ring 
; [denoted by (Be), in the ring model; (B,-), and (Be), in the ax- 
isymmetric model; brj, bgj and 0,- in the bisymmetric model]. 
Therefore, the 9 different rings may be analyzed separately. 

For any one of the three field models, consider a given ring / 
and suppose that this ring contains n, boxes. For every box j, we 
have derived an observational value of the average line-of-sight 
field (B||^, denoted by (^11)... together with its statistical uncer- 



tainty, denoted by (cr||),y. Besides, we can calculate a theoretical 
expression of the large-scale line-of-sight field (By), denoted by 
{B\\)jj, in terms of the free parameters of ring /. The best-fit val- 
ues of these parameters are obtained by minimizing 



X 



(^11)/; 



(13) 



The best fit of ring / can be considered consistent with the 
pulsar data if, on average over ring /, the theoretical best-fit {Bn)^. 

do not differ from the observational ^By^ , , by more than the asso- 
ciated uncertainties (crg)ij. In mathematical terms, this condition 
for consistency can be expressed asx^ ^ «/• However, when the 
number of data points, n,, is not much greater than the number 
of free parameters, v [v = 1 for the ring model; v - 2 for the 
axisymmetric model; v = 3 for the bisymmetric model], consis- 
tency with the data should be tested with the more exact criterion 



X < «/ ■ 



(14) 



where n, -vis the number of degrees of freedom, i.e., the number 
of data points that cannot automatically be placed on a curve 
with V adjustable parameters. Eq. ( fT4l i provide s a rule of thu mb 
for a reasonably good fit (see Sect. 15.1 in Press et al.lll992h . If 
» n, - V, the best-fit curve misses too many data points to be 
believable. 

It is possible to obtain a more rigorous (and, at the same 
time, more flexible) criterion for consistency. Suppose, for the 
sake of argument, that the model we are testing is correct. If the 
data points (By^ , of ring / follow a Gaussian distribution, has 
a chi-square distribution with n, - v degrees of freedom. One 
can then calculate the a priori probability, P\x^ > xlnt^^ 
the obtained for a particular set of data points exceeds some 
critical value Conversely, one can calculate the critical 
for which Plx^ > xlnJ equals some imposed probability Pq. For 
instance, if the model is correct, it is unlikely (only 10% chance) 
that;^'^ > xlntiPo - O-l)- Turning the statement around, if we 
find;i^^ > x^ XPo - 0.1), it is unlikely that the model is correct 
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Table 1. Critical values ofx^ for 3 probability levels 



Hi - V 


xlJPo = 0.05) 






1 


3.841 


2.706 


1.642 


2 


5.991 


4.605 


3.219 


3 


7.815 


6.251 


4.642 


4 


9.488 


7.779 


5.989 


5 


11.070 


9.236 


7.289 



- we will say that the model is inconsistent with the data. This 
reasoning directly leads to the following consistency condition: 



(15) 



Here, we will adopt Po - Q.l as our default value, but we will 
also discuss the results obtained for Pq - 0.05. In practice, the 
values of Plx^ > xlni] f or given w, - v and;^f^^.j^ are tabulated in 
various textbooks (e.g., lYamanelll964j) . The tables can also be 
used to determine ;t'crit*-^0-' given «, - v and Pq. For reference, 
the values of xl^itiPo) for n,- - v = 1, 2, 3, 4, 5 and for f o - 
0.05, 0.1, 0.2 are listed in Table[T] An important point emerging 
from Table[T]is thatxl^^^iPo = 0.1) > n, - v, so that Eq. ([151) with 
Pq - Q.l will always be easier to satisfy than Eq. (fT4b . 

In the next three subsections, we present the results obtained 
with the rule of thumb (Eq. (fT4l) ) and with the more rigorous 
consistency condition (Eq. ([TS])). for the three field models. 



4.1. Ring model 

In the ring model, (B, ) vanishes and (Bg) is constant along cir- 
cles. Hence, there are 9 free parameters: (Bg),, the large-scale 
azimuthal fields in the 9 rings / - 4, 4.5, 5 ... 8 Q Since all the 
free parameters are independent, the 9 rings can be treated sepa- 
rately. 

For every ring the large-scale line-of-sight field in any box 
j is simply the projection of (Be),- onto the line of sight (see 
Eq. © with (Br) = 0): 



(^ll>,7 = (Be)i cos Oij 



(16) 



where ajj is the angle between the azimuthal direction and the 
direction to the Sun at the midpoinQ of box j (see Fig. [1]). 
The best-fit value of (Be); is obtained by minimizing x^ (given 
by Eq. (fT3Tl). In terms of (fie).. - {^B\\j.. / cos atj, the obser- 
vational value of the average azimuthal field in box j, and 
(o"e),7 = (cr||),77 cosa,;, the associated uncertainty, the min- 
imization procedure turns out to be equivalent to taking an 
uncertainty-weighted average of the different (Bs).. along ring 
;: 



{Bg}i = 



f Hj 



(17) 



' As a reminder, ring / is centered on r; = ; kpc and extends between 
((■ - 0.5) kpc and (i + 0.5) kpc. 

^ As in Fig. [4] the midpoint of a box is defined as the point at middle 
radius, r = r,, and middle longitude, / = + /max)/2. 



The values of (Be)., and their uncertainties (crg)ij in the n, 

boxes j of the 9 rings / are plotted in Fig.|5] at the Galactic az- 
imuthal angles of the box midpoints, 6,7. For comparison, the 
best-fit values of (Be), in the 9 rings are indicated by horizontal 
lines spanning the entire azimuthal range. 

In only one ring (; - 4.5) does the best-fit value of (Be),- 
satisfy the rule-of-thumb consistency condition, x^ ^ ^ 1 
(Eq. (fT4l i with v = 1). For this ring, we compute the consistency 
range of (Be),, which contains all the values of (Bg),- for which 
X^ < rtj - I. The best-fit value of (Be), and its consistency range 
in the sole "good-fit" ring are plotted against r,, in the upper 
panel of Fig. |6] For the other 8 rings, the (inconsistent) best-fit 
values of (Bg), are plotted with crosses. Clearly, these 8 rings do 
not admit any ring magnetic field consistent with the data. As an 
immediate consequence, the ring model must be rejected. 

In order to gain some feel for how far the ring model is from 
being able to reproduce the pulsar data, let us, in thought, extend 
the error bars of all the observational (Bh^ . by a factor of 2 and 

look into the impact of this extension on our results. With twice 
the original error bars, the parameter would be smaller by a 
factor of 4, so that, in terms of the original the rule-of-thumb 
consistency condition would become < 4(n,- - 1). As it turns 
out, this less stringent consistency condition would be fulfilled 
in 5 rings (/ = 4, 4.5, 5.5, 6, 7) out of 9. Thus, with twice the 
original error bars, the ring model would remain unacceptable. 

If we now resort to the more rigorous consistency condition, 
X^ ^ xlnt(0-l) (Eq. (O with Pq - 0.1), to test the ring model, 
we find that 3 rings (/ - 4, 4.5, 6) have their best-fit (Be), con- 
sistent with the data; their consistency ranges are drawn in the 
lower panel of Fig. |6] For the other 6 rings, the (inconsistent) 
best-fit (Be),- are again plotted with crosses. With Po = 0.05, 
4 rings (/ = 4, 4.5, 6, 7) would be deemed consistent with the 
data, but the other 5 rings would still fail the consistency test. 

These results confirm our conclusion that the ring model 
must be rejected. 

4.2. Axisymmetric model 

In the axisymmetric model, (B^) and (Bg) are both constant 
along circles. Hence, there are 18 free parameters: (B^), and 
(Be),, the large-scale radial and azimuthal fields in the 9 rings 
/ = 4, 4.5, 5. ..8. 

For every ring /, the large-scale line-of-sight field in any box 
j can be written as a linear combination of the two parameters 
(B,),- and (Be),-: 



(B||).. = (Br),- sinQ',7 + {Bg)i cos aij 



(18) 



(see Eq. (|5]l). Again the best-fit values of (Br), and (Be), are ob- 
tained by minimizing (Eq. (fT3]l). 

Here, we find that 1 ring (/ = 4.5) has its best fit consistent 
with the data, according to the rule-of-thumb consistency condi- 
tion, x^ < n,- - 2 (Eq. ( fT4l l with v = 2). Its consistency domain 
in the parameter plane ((Br), , (Be),) is the area delimited by the 
ellipse x^ - «/ - 2 (grey contour line in the relevant panel of 
Fig.|7J. As none of the other 8 rings can be properly fit with an 
axisymmetric magnetic field, the axisymmetric model must be 
rejected. 

With twice the original error bars on the observational 

the rule-of-thumb consistency condition would become x^ ^ 
4(n, - 2) (in terms of the original ;^'^). This less stringent consis- 
tency condition would be satisfied in 5 rings (/ - 4, 4.5, 6, 7, 8), 



H. MEN et al.: Observational constraints on models for the Galactic magnetic field 



7 



ring 4 



ring 4.5 



ring 5 



ring 5.5 



ring 6 



ring 6.5 



ring 7 



ring 8 



a. 



a. 



la I 





■ 1 1 1 ■ 1 _ 

i.L'.-'.u.i. . - 







a. 



o 



a. 



a. 



loT 4 



i 



-.■-,■4 





r ' _ 







1 . ■ . 1 


■ 1 1 1 ■ 1 







7.5 










' _ 





CN 

4 
1 






H 1 1 1 ] 

1 . . " 














8 










1 1 1 1 _ 

H ^ ^ — : 




-4-2 C 






1 . . " 










-50 

e.. (deg) 


50 



Fig. 5. Observational values of the average azimuthal fields, (^Bg^ . ,, veriMi Galactic azimuthal angles, Gij, in all the boxes /' of the 9 

successive rings /, for the ring model. The uncertainties (crg),-^ in the field values are plotted as standard error bars. For each ring 
the best-fit value of the model parameter (Be), is indicated by the horizontal dot-dashed line. 



but still not in the other 4 rings. Therefore, the axisymmetric 
model would remain unacceptable. 

According to the more rigorous consistency condition, < 
xlj(^-l) (Eq. dEl) with Po = 0.1), 5 rings (/ = 4, 4.5, 6, 7, 8) 
have their best fits consistent with the data. Their consistency 
domains in the parameter planes ((Br); , (Be)/) are the elliptical 
areas enclosed by the curves - ^cnt*-^-^-* (black contour lines 
in the relevant panels of Fig. |7]l. For the other 4 rings, the (in- 
consistent) best-fit pairs ((B, ),- , (Bg),) are indicated with crosses. 
Relaxing the probability level to Pq = 0.05 would not raise the 
number of acceptable rings above 5. 

From all the above, we conclude that the axisymmetric 
model must be rejected. 

4.3. Bisymmetric model 

In the bisymmetric model, (B^) and (Bg) vary sinusoidally along 
circles in the manner described by Eqs. (fTOl i and ( fTTT i. Hence, 
there are 27 free parameters: brj, bgj and 0,-, the maximum am- 
plitudes and the azimuthal phases in the 9 rings / = 4, 4.5, 5 ... 8. 



For every ring /, the large-scale line-of-sight field in any box 
j follows from Eq. dD together with Eqs. (fTOli- (fTTl i: 

(B||),.^. = brj sin(% - (pi) sinaij + bgj sm(6ij - (pi) cosa,y , (19) 

where the angles Oij and aij (see Fig. [T]l refer to the midpoint 
of box /'. Similarly to the previous models, the best-fit values 
of the three parameters b, j, bgj and are obtained through a 
minimization of;^^ (Eq. (fTSl)). 

According to the rule-of-thumb consistency condition, ^ 
«, - 3 (Eq. ( fT4l i with v - 3), 2 rings (/ - 4.5, 5) have their 
best fits consistent with the data. Their consistency domains in 
the parameter spaces (brj, bgj, (pi) are the volumes bounded by 
the surfaces = n,- - 3. Displayed in Fig.|8]are the projections 
of these consistency domains on the parameter planes {brj, bgj) 
(grey contour lines). Since the other 7 rings fail the consistency 
test, the bisymmetric model must be rejected. 

With twice the original error bars on the observational 
{B\\)..7 the rule-of-thumb consistency condition would become 

< 4(n,- - 3), which would be satisfied in 7 rings (/ — 
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Fig. 6. Best-fit values of the model parameters, (Be),, versus ring middle radii, r,, for the ring model. The acceptable values (those 
consistent with the pulsar data) are shown with their consistency ranges, while the unacceptable values are indicated with crosses. 
The results obtained with the rule-of-thumb consistency condition, < n, - 1, are plotted in the upper panel. Those obtained with 
the more rigorous consistency condition, ^ ,(0. 1), are plotted in the lower panel. 



4, 4.5, 5, 5.5, 6, 7, 8) out of 9. The bisymmetric model would 
then be nearly acceptable. 

According to the more rigorous consistency condition, ^ 
xlj(^-l) (Eq. O with Po = 0.1), 5 rings (/ = 4.5, 5, 6, 7, 8) 
have their best fits consistent with the data. Their consistency do- 
mains, bounded by the surfaces = ^crit^'^-^-'' shown 
in projection on the parameter planes (brj,bgj) in Fig. [8] (black 
contour lines). For the other 4 rings, the (inconsistent) best-fit 
pairs (b,- i, bgj) are indicated with crosses. Relaxing the proba- 
bility level to Po = 0.05 would raise the number of acceptable 
rings to 7 (/ = 4, 4.5, 5, 5.5, 6, 7, 8), which would render the 
bisymmetric model nearly globally acceptable. 

Altogether, the bisymmetric model must be rejected, though 
its rejection is slightly less severe than for the axisymmetric 
model. 

5. Summary and conclusions 

In this paper, we examined the three most common theoretical 
models for the large-scale magnetic field in the Galactic disk 
and confronted each of these models with the pulsar data. For 
each model, we derived the best-fit parameters, through min- 
imization, in the 9 Galactocentric rings defined in Fig.|2l and we 
delineated the parameter domains around the best fits (referred 
to as the consistency domains) wherein the predicted fields are 
consistent with the pulsar data. 

Compared to existing studies of the kind, we did not attempt 
to settle the long-standing (and possibly ill-posed) question of 
whether the Galactic magnetic field is axisymmetric or bisym- 
metric. Our sole purpose was to determine whether each of the 
three basic field models, taken separately, is compatible with the 
available pulsar data or not. In this regard, we note that many 
previous studies did find a preference for one of the field mod- 
els, but omitted to put their preferred model through the crucial 



"goodness-of-fit" test, which checks whether the model can in- 
deed reproduce the data within the error bars. 

Here, we tested the three field models on the basis of two 
different criteria: first, a standard rule of thumb for a reason- 
ably good fit (Eq. (fT4b ). and second, a more rigorous consistency 
condition for a chi-square distribution of (Eq. ( fTSl ). with the 
imposed probability set to Pq - 0.1). These two criteria were 
successively applied to all the rings separately, such that the best 
fit of ring / was deemed consistent with the pulsar data if the as- 
sociated value oix^,X^in' was less than n,- - v (first criterion) or 
less than xl^-i^iPo = 0.1) (second criterion). A model could then 
be considered globally acceptable if all the rings had their best 
fits consistent with the data. 

The results obtained for the three field models, with both cri- 
teria, are summarized in Table |2] All the rings are listed with 
their labels, / (see footnote [1]), and their numbers of boxes (or 
numbers of data points), n,, from which it is straightforward to 
deduce the numbers of degrees of freedom, n, - v (for the first 
criterion), and the critical ;\f^j.jj(Po = 0.1) (for the second crite- 
rion; see Table [1]). Also given for all the rings are the minimum 
values ofx^, ^Ln' '■^^ values associated with the best fits, as 
well as the results of both consistency tests (satisfaction of a test 
is indicated with an asterisk), for the three field models. 

We found that none of the three field models is acceptable, 
in the sense that none of them can be brought into full agree- 
ment with the pulsar data. According to the standard rule-of- 
thumb consistency condition (Eq. (fT4b ). all three models must 
be strongly rejected, as the ring and axisymmetric models fail to 
provide a good fit (consistent with the data) in all the rings save 
one, while the bisymmetric model fails in all the rings save two. 
If the error bars of the observational line-of-sight fields were en- 
larged by a factor of 2, the bisymmetric model would not be too 
far from acceptable (7 good-fit rings out of 9), while the ring and 
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Table 2. Summary of the results obtained for the three field 
models " 



Ring model (v = 1) 



Axisymmetric model (y = 2) 



Bisymmetric model (v = 3) 



i rii 


mm 


4 


4 


5.03 


4.5 


6 


3.17 


5 


5 


25.37 


5.5 


5 


14.52 


6 


6 


8.84 


6.5 


5 


20.29 


7 


6 


9.57 


7.5 


4 


19.31 


8 


4 


52.53 



^ mm ' 



y2. < 1/2.(0.1) 
•^min — /vcnt^ 



■ <ni-y . <y^ .,(0.1) 

mm ' mm crit ^ 



X - <rii 

mm ' 



•^min — -^crit^^" ^-^ 



2.91 

3.11 
18.65 
14.43 

5.06 
17.93 

5.81 
19.30 

3.04 



3.08 
1.75 
1.79 
4.95 
5.48 
13.55 
3.37 
7.44 
0.82 



" When a value of x~ ■ satisfies the first or second consistency condition, an asterisk is plotted in the corresponding column. 



axisymmetric models would remain truly unacceptable (5 good- 
fit rings out of 9). The conclusions reached with the more rig- 
orous consistency condition (Eq. (fTSl l) are intermediate between 
the two situations described above; with the imposed probability 



set to Po -0.1 (Po - 0.05), 3 (4) rings can be properly fit with 
a ring magnetic field, 5 (5) with an axisymmetric field and 5 (7) 
with a bisymmetric field. 
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The quantitative differences between both criteria are eas- 
ily understood. A comparison between Eqs. ( fT4l i and (fTsT i 
immediately shows that the rule-of-thumb consistency condi- 
tion (Eq. ( fT4l l) corresponds to a probability level Pi such that 



or equivalently. Pi - P\x > n-i - v\. With 



twice the original error bars on the observational line-of-sight 
fields, the rule-of-thumb consistency condition would become 
< ^(n-i - v), corresponding to a probability level P2 such that 
xiitiPi) = 4(n,- - v), or equivalently, P2 = Plx^ > Mrii - v)]. The 
values of Pi and P2 as functions of n, -v are tabulated in Table|3] 
Clearly, the rule-of-thumb consistency condition with the origi- 
nal error bars implies high probability levels (Pi ~ 30% - 40%), 
which make it overly diflicult to satisfy; if the model under test- 
ing is correct, there is nonetheless a ~ 30% -40% chance that;i'^ 
exceeds n, - v and that the model will be rejected. On the other 
hand, with twice the original error bars, the probability levels 
drop very low (P2 ~ 0.1% - 4%) and the consistency condition 
becomes too easily satisfied; the risk is then to accept a model 
that is in fact incorrect. 

The bottom line is that the standard rule of thumb is way too 
stringent, while the rule of thumb with twice the original error 



Table 3. Probability levels for the rule-of-thumb consistency 
condition 

n, -V I Fi=Plx^> -y] P2 = P\X^ > 4(«, - y)] 



0.3173 
0.3679 
0.3916 
0.4060 
0.4159 



0.0455 
0.0183 
0.0074 
0.0030 
0.0012 



bars is way too easy. In contrast, the more rigorous criterion with 
a probability level set to Pq = 0.1, intermediate between Pi and 
P2, provides a reasonable trade-off. The latter criterion is also 
more trustworthy, insofar as all the rings for all the models are 
tested with the same probability level. 

The results of the present study suggest that the true large- 
scale magnetic field in our Galaxy has a more complex config- 
uration than a strictly axisymmetric or bisymmetric field. One 
possibility is that it consists of the superposition of axisymmet- 
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lie, bisymmetric and probably higher-order azimuthal modes. 
Such combinations of modes have been observed in a number of 
external galaxies ( e.g., Beck et al.,1996i ; ,Berkhuiisen et al..l997t 
iRohde et al.lll999HBeckll2007h . 
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